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CHARACTERIZING THE SCIENTIFIC POTENTIAL OF 
SATELLITE SENSORS 
GSFC ATR - Dr. J. Barker 
SAR Ta^ Leader - Dr. Y Lee 


Taric Objective : 

The objective of this task is to provide analytical and programming support to characterize 
the potential of the Landsat Thematic Mapper (TM) digital imagery for scientilic 
investigations in the earth sciences and in terrestrial phykcs. Secondary objectives of this 
task include providing technical support to define lower atmospheric axxl terrestrial surface 
experiments for the space station and providing technical sunxnt 
Sensor (ROS) stucty scientist for advanced studies in remote sensing. 


Work Performed : 

Tne following work was performed in the areas indicated. 

100 TM Radiometric Calibration 

Eleven TM radiometric calibration programs developed last montlhfa fefe bfe tjl 4[^ed 
and evaluated. Several software problems have been found and^»h^. The 
correction algorithm used for coherent noise correction has been modified. 
Unsuc::essful correction of bright target saturation in the background region h<»s 
been examined. 

Three software errors have been found and solved as follows. 


> the Research Optical 



• Integer Overflow - Because the VAX integer*2 format used in the calibration 
has a sign bit. the masdmum value cannot be 65.535. The problem was solved 
by using 14-bit calibration (maximum value 16.383) instead of 16-bit 
calibration. The original 8-bit TM data then is expanded approximately 64 
times. 


• Divide by Zero - During the calibration of masks for scan-correlated shifts, a 
threshold value for separating lower and higher states was calculated using ^e 
data from both forward and reverse scans. This technique of threi^olding may 
deflect either all forward scans or all reverse scans as being in either a higher 
or lower state. The average value of each state is calculated by summing the 
values for each state and dividing by the state population. When all forward or 
reverse scans are in one state, the opposite state will have a zero population 
and the divide-by-zero problem will occur. This problem has been solved by 
defining separate threshold values for forvard arul reverse scans. 
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• Nonexist File Group - This error was caxised by mistakenly defining the image 
I/O file name in several programs. It has been corrected. 

Since most of the output files from the 11 calibration programs are in 14-bit 
format, it is necessary to convert tlxe files into t^e format for display on the 
Interactive Analysis Terminal (lAT). A program. IBCON.PDF. has been implemented 
for this purpose. It performs a linear stretch between the minimum and maximum 
valiies to convert 14-bit image data into 8-bit data. 

The algorithm for coherent noise correction has been examined. The results from 
the scan-by-scan correction procedure utlU^ng single-band resequenced data did 
not ^ow a consistent characteristic through scans in the band 1. channel 9. noise 
image. The line-by-line correction procedure suggested by the ATR has been 
implemented. 

A problem in applying a bri^.t target saturation correction in the background region 
is under examination. It may be due to error in the boundary mask or the recovery 
time constants. 

200 Image Processing on LAS/VAX 

The following raw. calibrated, aixl corrected image groiq^s have been created and 
stored on the Barker2 disk. These image groups may be deleted and recreated after 
applying the BTS correction in the background region and examining the results. 

• CSF4 - The raw image group includes images of bands 1 through 4 of a 

512 X 12 sbbscene at (SP. SL) » (2101. 2337) of scene 40392-18152 
over San Francisco. 

• CBTS - BTS corrected Images of the CSF4 group. 

• CSCS - BTS+SCS corrected images of the CSF4 groiQ). 

• CCN - BTS+SCS+CN corrected images of the CSF4 group. 

• CCAL - 8TS+SCS+CN+TRAPP (IC+H) calibrated images of the CSF4 

group. 

230 Image Processing on IDIMS/HP 

Black-and-white pixel print files were created for the following San Francisco 
subscenes (Scene ID 40392-18152): 
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Name 

Band 

SP. SL. NP. NL 

Processina Information 

SFB.Bl 

1 

6629. 1. 294. 384 

Background data without 

corrections 

SFB.CNl.L 

1 

6629. 1. 294. 384 

Background data with coherent 
noise correction applied on a 
line- by-line basis 

SFB.CNl.S 

1 

6629. 1. 294. 384 

Background data with coherent 
noise correction applied on a 
scan-by-scan basis 

NOISEl.L 

1 

6629. 1. 294. 384 

Noise image from line-by-line 
method 

NOISEl.S 

1 

6629. 1. 294. 384 

Noise image from scan-by-scan 
method 

NOISEl.16.9 

1 

6629. 1. 294. 384 

Noise image for 32 scans of 
channels 9 and 16 from 
sc«:i-by-scan and line-by-line 
corrections 

SFM.Bl 

1 

1825. 1889. 294, 384 

No corrections applied 

SFM.Bl’ 

2 

182.5. 1889. 294, 384 

No corrections applied 

SFM.B3 

3 

1825, 1889. 294, 384 

No corrections applied 

SFM.B4 

4 

1825. 1889, 294. 384 

No corrections applied 

SFM.BTl.T 

1 

1825, 1889, 294. 384 

BTS correction using truncation 

SFM.BT2.T 

2 

1825. 1889. 294. 384 

BTS correction using truncation 

Black-and-white negatives were made and sent to the i^otographic laboratory to 
have view graphs and five prints made of each. 

These products were reviewed by the ATR and found to have low contrast over 
water and high contrast over land. 
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300 Production Support 

Becasue of a lack of source code of dump utilities for the fUe on the CALFILE tape, 
it is difficult to generate data from the tape (for example, the housekeeping data) in 
the digital form. The CALFILE processing will continue as before, producing a 
photocopied output book of the following files: TAG. (Quality Assurance Data File 
(QAD), Housekeeping Data File (HSD), Ephemeris Data File (ESD), Systematic 
Correction Data File (SCD), Control P(^t Data File (CPD), Long-Term Parameter 
File (LT-PAR), Short-Term Parameter File (ST-PAR), Histogram File (HISTPLT). 
Scan Characteristics Record (SCH), and Reduced Calibration File (CAL). A 
specified data file may be created independently if requested by the ATR by 
manually copying sections from the photocopied output book or editing the output 
file stored on the VAX system. 

The incorrect telemetry data on the 45 BRU tapes received from SBRC have been 
confirmed with B. Cooley at SBRC. The correct data have been recovered from 
A\ig. 30, 1983, BARVID printouts. SBRC will examine problems of BURST data arid 
data from bands 5, 6, and 7. No further action will be taken on processing the 45 
tapes until data verification is received from SBRC. 

Six CALFILE tapes and one CALDUMP tape were processed for the following scene 
IDs: 


ID 


4060815463 

5001415460 

4039218152 

5012917075 

5014216365 


Location 

Birminfhani 
Birmingham 
San Francisco 
White Sands 
Unknown 


Four blue books were generated for each scene. 

Two chosen SBRC tapes (#1038. 5 window, aM #10/1, 8 window) have been run 
successfully. 

Two BURST files were retrieved from the SBRC tapes (#0138, #1039) and have been 
processed with TRAPP. The full scan line was plotted for each file. 

All the BRU data files from the 45 SBRC tapes have been copied and merged onto 
eight BRU multifile tapes - #309, #310, #311, #312, #313, #314, #315, #316. 

Twenty-one calibrated and corrected subimages of San Francisco were created and 
transferred onto tape for producing pseudocolor pixel prints. 
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400 Software Development 

Linear regression software is being dc/eloped to correlate the start of scan time 
with the shutter background (Section 410) and to correlate background data between 
channels within a band (Section 420) and between channels in different bands 
(Section 430). Since there should be no correlations in the background data in the 
absence of systematic error s , the outputs from these programs are designed to 1) 
investigate the origin of systematic radiometric variability, such as coherent noise 
and scan-correlated shift, and 2) serve as figures-of-merit to evaluate the 
effectiveness of correction for these systematic errors. 

410 Linear Regression Coefficient Tables, Part I 

This software performs the following optional regressions of shutter background 
digital number (DN) versus: 

• Start of scan in seconds (PO) 

• Difference between adjacent scans of the start of scan in seconds (DPO) 

• Difference of DPO from the average val\ie of DPO in seconds (DDPO) 

The output from this software consists of four pages, as follows: 

The first page is a list of pairs of avenged shutter background (BOBS) and either PO, 
DPO, or DDPO by scan line for the selected band and channel (see Appendix 1, 
results of Landsat-4 BRU tape GSFC 4 :?). 

The second page is a list of backgound residuals (BCALC-BOBS) calculated from 
linear regre^on of BOBS versus PO, DPO, or DDPO by scan line for the selected 
band and cV^annel (see Appendix 2, resu^ ts of Landsat-4 BRU tape GSFC #37). 

Pages i and 4 contain eight tables of derived values for all 100 charjiels (see 
Appeniiixes 3 and 4, results of Landsat-4 BRU tape GSFC #37). 

• AO(DN). intercept of linear regression 

• SDA0(DN). standard deviation of the intercept 

• Al(DN/sec), slope of linear regressions 

• SDAl(DNZsec), standard deviation of the slope 

• R, correlation coefficient 

• Randomness of background residuals 

• SB(DN). standard error of background 

• 100*(1-R2), unaccounted variances in percent 

All the above outputs are in floating point format as requested by the ATR, and 
have become part of TRAPP’s output. A brief description of this software 
development can be found in the previous monthly report; for the details, contact 
P. Lee, SAR. 


5/8 


Ta^ Assignment 132 
December 1984 


420 Linear Regression Coefficient Tables. Part n 

This part of the analyses study channels correlations within a band by regressing the 
^tter background data for channel I versus the ^tter back^xnmd data for 
channel J. Five different optional output can be chosen, as follows: 

• All scans versus all other scans 

• Even scans versus even scans 

• Odd scans versus odd scans 

• Even scans versus odd scans 

• All of the above 

Therefore, a maximum of 28 pages may be generated for the 7-band data whei the 
final option above is chosen. Each page of output contains the following three 
tables: 

• Correlation coefficient, R, and the unaccounted variance. 100 (1-R^2) 

• The intercept and standard deviations of the intercept 

• The slope and standard deviation of the slope 

The following two new routines have been implemented and tested: 

• Subroutine RCORR02 - This routine generates the output tables for the 
within-band regressions. 

• Subroutine RCFORM - This routine formats all output variables from the 
within-band regression 

Minor changes in routines RCINJBRU, RCUSERIN, and RCNMAIN to link the above 
new routine to the TRAPP subroutine library m^dee these new tables part of the 
standard TRAPP output. 

430 Linear Regression Coefficient Tables. Part III 

This part of the analyses study channels correlations between bands by regressing 
the Gutter background data for reference channels of each band. The software 
outputs three pages and nine tables. The new routines imder development are as 
follows: 

• Subroutine RCCORR03 - This routine generates output tables for the 
between-band regressions. 

• Subroutine RCOUT - This routine formats the output variables. 

• Subroutine RCIN - This routine allows the user to prompt the input and ^ows 
it on a screen for double checking. 
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700 Image Science Studies 

The contouring program on the I2S was tested by the ATR. It was found that a 
512 X 512 image takes 30 minutes to process and uses 80 percent of the VAX 11/780 
CPU. 

TM images of Jiddah. Saudi Arabia (Scene ID 5009707171. path 170, row 45. 
June 6, 1984), and northern Chile (Scene ID 5012914060, path 1, row 77, 
July 8. 1984) were ^^'wn to the ATR on the Land Analysis System (LAS). The 
Jiddah scene has random speckle noise in all bands, and the Chile scene has response 
differences between back and fore scans in band 7. 

The ATR and SAR discussed the following questions concerning the image science 
studies: 

1. What are the scientific objectives of (\oing corrections to TM data %’ersus the 
applications? 

2. What type of parametric error analysis can SAR perform to quantitatively 
assess the improvement in TM image data quality after doing radiometric 
corrections? 

3. How can the impact of these improvements be demonstrated in science 
studies? 

All results for the January TM workshop for science studies will be demonstrated 
using the San Francisco scene. 


Problem Areas : 

230 Image Processing on IDIMS/HP 

No image processing was performed on the HP-3000 during December 10-21 because 
the Laboratory for Terrestrial Physics Computing Facility is under reconstruction. 

Schedule Conformance : 

Work is proceeding as planned. 

Work Planned for Next Month: 


100 TM Radiometric Calibration 

Task personnel will reexamine parameters for bright target saturation. 
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200 Image Processing on LAS/VAX 

Calibrated and corrected image and background files and the difference image files 
will be created. 

230 Image Processing on IDIMS/HP 

Pixel prints will be created of image and background data files for various types of 
calibrations. 

300 Production Support 

SAR personnel will i^xlate apparent gain versus time plot for both Landsat-4 and -5 
TM. 

SAR personnel will continue to process new CALDUMP and CALFILE tapes (6 sets). 

400 Software Development 

Major emphasis wi i be placed on analyzing and implementing the different background 
linear regressions. A task member will continue to implement and test between-band 
linear regicr<tion coefficients and make any format changes to the previously 
developed software upon the ATR’s request. 

700 Image Science Studies 

Work will start on assessing the impact upon information extraction of doing various 
types of radiometric calibrations to the data. 


Deliverables Submitted : 

Graphics: 12 black-and-white negative pixel print films of raw and corrected TM images 

and background data 
Originator: J. Wang and W. Hallada 


Computer Utilization : 

The estimated computer time used this month is as follows: 


Minutes 

1501 (wall clock) 
200 (wall clock) 


Computer 

HP-3000 (IDIMS) 
VAX 11/780 (LAS) 
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SUBTASK »00 

FINANCIAL DATA AND LABOR-HOUR UTILIZATION 


Financial Data 

Cumulative total to date: 
Projected total to complete: 
CTR estimated total: 


$ 33.253 
$ 90. 923 
t 124.1 76 


Labor-Hour Utilization 


Hours this Month 


Cateaorv 

On-Site 

Off-Site 

PM 


3.6 

TAM 


9.0 

LPA 

19.5 

18.0 

SAP 

40.0 


AP 

eo.o 


JAP 

160.? 


SP 

158.0 


Support 


30.7 


Cumulative Hours 
On-Site Off-Site 

11.8 

32.0 

64.0 68.5 

137.0 
276.0^ 
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